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We show how microstructure can arise in first-order ferroelastic structural transitions, in two and 
three spatial dimensions, through a local meanfield approximation of their pseudospin hamiltoni- 
ans, that include anisotropic elastic interactions. Such transitions have symmetry-selected physical 
strains as their Aop-component order parameters, with Landau free energies that have a single zero- 
strain 'austenite' minimum at high temperatures, and spontaneous-strain 'martensite' minima of Ny 
structural variants at low temperatures. The total free energy also has gradient terms, and powerlaw 
anisotropic effective interactions, induced by 'no-dislocation' St Venant compatibility constraints. 
In a reduced description, the strains at Landau minima induce temperature-dependent, clock-like 
Zjv v +i hamiltonians, with TVop-component strain-pseudospin vectors S pointing to Ny + 1 discrete 
values (including zero). We study elastic texturing in five such first-order structural transitions 
through a local meanfield approximation of their pseudospin hamiltonians, that include the power- 
law interactions. As a prototype, we consider the two- variant square/rectangle transition, with a 
one-component, pseudospin taking Ny + 1 = 3 values of S — 0, ±1, as in a generalized Blume-Capel 
model. We then consider transitions with two-component (Nop — 2) pseudospins: the equilateral 
to centred- rectangle (Ny = 3); the square to oblique polygon (Ny = 4); the triangle to oblique 
(Ny = 6) transitions; and finally the 3D cubic to tetragonal transition (Ny = 3). The local mean- 
field solutions in 2D and 3D yield oriented domain-walls patterns as from continuous-variable strain 
dynamics, showing the discrete-variable models capture the essential ferroelastic texturings. Other 
related hamiltonians illustrate that structural-transitions in materials science can be the source of 
interesting spin models in statistical mechanics. 

PACS numbers: 81.30.Kf, 64.70.Nd, 05.50.+q, 75.10.Hk 



I. INTRODUCTION 



Ferroelastic crystals undergo diffusionless structural 
transitions that are first order, and on cooling show a 
reduction in symmetry to two or more spontaneously- 
strained states (or Variants') which can be transformed 
between one another by stress^. These transitions are 
often studied through minimizing the Landau free ener- 
gies 3 - in terms of appropriate continuous variables, such 
as displacements, phase fields, or strains^—. Although 
homogeneous, single- variant martensite states are the 
global minimum, clastic heterogeneities or mctastablc 
domain-wall patterns are experimentally found^, that are 
locked in to preferred crystallographic directions^. This 
orientation arises through a balance between the Lan- 
dau energies nonlinear in the order parameter strain, 
the short-range gradient costs or Ginzburg energies, and 
the effectively long-range elastic energies^ or powerlaw 
anisotropic interactions, that orient the domain walls 
J±±-1LL. The powerlaw interactions result from enforcing St 
Venant 'compatibility' constraint ;) 1 ^ 12 between the strain 
components so that the displacements are continuous, 
with no dislocations generated on cooling. Continuous- 
variable models have been used to study microstruc- 
tures under various conditions, including strain-rate de- 
pendencei 3 -; and the effects of finite-size on martensitic 
growth in an austenitic matrix—. 



Models in terms of discrete structure-variables or 
'pseudospins' have also been used to study these fer- 
roelastic transition a 15 ' 16 . (Similar in spirit to discrete 
strains, analytic 'minimizing sequences' consider tent- 
like displacement profiles or flat strain-variants, on ei- 
ther sides of domain walls 2 ..) Recently, model pseu- 
dospin hamiltonians induced by the scaled free energies 
for several specific transitions in 2D and 3D have been 
proposed^ 7 -. The model hamiltonian is simply the to- 
tal scaled^ free energy evaluated at the Landau minima 
in the order parameters (OP). The pseudospins are 'ar- 
rows' in A^op-dimcnsional order-parameter space, point- 
ing to the Ny variant minima, and to the zero-strain 
turning point. The hamiltonian includes a temperature- 
dependent on-site term quadratic in the pseudospins from 
the Landau term, a nearest-neighbor ferromagnetic inter- 
action between pseudospins from the Ginzburg term, and 
a pseudospin powerlaw interaction from the St Venant 
term. The pseudo-spin models are like Zjv v +i clock mod- 
els^ generalized to include a spin-zero state, and may be 
termed 'clock-zero' models^ 7 -. A three-state spin-1 type 
model for the transition of square to rectangle unit cells 
(with Nop = 1,AV = 2) has found glass-like behavior 
on slow cooling using a local meanfield approximation^ 7 -. 

In this work we consider pseudo-spin models in 
the local meanfield approximation under temperature 
quenches, for five structural transitions : four in two 
spatial dimensions 2 ^, and one in three spatial dimcn- 
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sionsi. Apart from the single order parameter (Nop = 1) 
square/rectangle case that is first studied as a simple pro- 
totype, the other four transitions all have two-component 
(Nop = 2) pseudospins. The transitions are i) the square 
to rectangle (2D version of the tetragonal/orthorhombic 
transition such as in YBCO); ii) the square to oblique 
polygon; iii) the triangle to centred-rectangle (2D ver- 
sion of hexagonal to orthorhombic transition such as in 
lead orthovanadate^i) ; iv) the triangle to oblique; and 
v) in 3D, the cubic to tetragonal transition (as in FePd). 
For these five transitions, the (nonzero) pseudospin ar- 
rows point respectively, to the two ends of a line, and 
to the corners of a square, a triangle, a hexagon, and a 
triangle, with the number of pseudospin variant states 
thus being respectively Ny = 2,4,3,6 and 3. We show 
that these discrete-variable models, despite their sim- 
plicity, produce local-meanfield microstructure in onc- 
and two-component strain pseudospins in agreement with 
continuous-variable strain simulations*—, that can be 
computationally more intensive. We thus find paral- 
lel twins for the square/rectangle and cubic/tetragonal 
transitions; nested stars for the equilateral/isosceles tri- 
angle transition; and tilted oblique domains for the 
square/oblique, and triangle/oblique transitions. 

The generalized clock model of strain pseudospins is 
a statistical mechanics description of the ferroelastic 
transitions in materials science. It conceptually links 
long-lived, metastable martensitic twins (even without 
quenched disorder) to Potts-model and clock-model de- 
scriptions of glasses^, and may be relevant to recent 
quenched-disorder strain glass behavior in martensitic al- 
loysSi. 

The plan of the paper is as follows. In Section |TT] 
we outline the derivation^ of the pseudospin hamilto- 
nians, and of compatibility potentials for the four 2D 
transitions. Our meanfield microstructure results are 
in Section IIIII where we first consider the two-variant 
square/rectangle case as a prototype, its response to ex- 
ternal stress, and its relation to the spin-1 Blume-Capcl 
model^i. We then consider local meanfield microstruc- 
ture for the three-variant triangle to centred-rectangle 
transition; the four- variant square/oblique transition; 
and the six- variant triangle/oblique transition. Turn- 
ing to 3D, Section ITVl considers the local meanfield mi- 
crostructure for the three- variant cubic/tetragonal tran- 
sition, with its compatibility potential stated in the Ap- 
pendix. In Scction[V]wc mention other related spin mod- 
els of interest in statistical mechanics. The final Section 
IVII has a summary and conclusion. 



II. PSEUDOSPIN HAMILTONIANS IN TWO 
SPATIAL DIMENSIONS 

The free energy functionals describing ferroelastic 
structural transitions can be written in terms of the phys- 
ical strains, that are symmetry-specific linear combina- 
tions of the Cartesian strain-tensor components. The 



Landau terms are invariant polynomials of the Nop 
order parameter strains, and have Ny minima. The 
free energies have many material-dependent elastic coef- 
ficients, that are not always known, or are fitted to exper- 
iment only for specific materials. However, the sponta- 
neous order-parameter strain magnitude at the first-order 
transition temperature is a small parameter. Following 
Barsch and Krumhansli^ a scaling procedure has been 
applied^ to four 3D transitions and five 2D transitions 
to obtain scaled Landau free energies that (to leading 
order in the small parameter) show universality at their 
minima, where any internal elastic constants are scaled 
out, and material dependence is only in an overall elastic- 
energy prefactor. 'Geometric nonlincarities' are higher 
order in the spontaneous strain and are neglected, as a 
perturbative first approximation. Then different mate- 
rials with the same transition, fall into the same 'quasi- 
universality' class, with common behaviour at the scaled 
minima, that lie at the corners and centres of the same 
'polyhedron' in Nop dimensional order parameter space. 
This is useful in strain- variable dynamics. It also im- 
mediately suggests a reduced description of fcrroclastics, 
in terms of discrete-strain statistical variables or vector 
'pseudospins', directed to these minima. 

A specific reduction procedure was proposed^ to ob- 
tain pseudospin hamiltonians by evaluating scaled free 
energies evaluated at their Landau minima. The ba- 
sic idea is quite simple, i) Scale the total free energy 
to dimensionless form, including the specifically calcu- 
lated compatibility-induced powerlaw interaction term, 
and the gradient term. Write the Landau free energy in 
polar coordinates in OP space, with the austenite min- 
imum at the origin, and Ny martensite minima located 
on a circle in Ny discrete angular directions, ii) Set 
the radial OP-magnitudc to its common temperature- 
dependent Landau-minimum value, and replace the OP- 
minima directional angles by discrete vectors pointing to 
these Ny + 1 minima on the circle, and at the centre, 
iii) The total free energy evaluated at minima is then 
the model Hamiltonian for the vector pseudospins, that 
have Nop spin-components, and Ny + 1 values. The 
remaining model coefficients are then not just arbitrary, 
but are related through the parent free energy, to the 
scaled temperature, to the scaled energy cost of an elas- 
tic domain-wall segment, and to the scaled bulk stiffness. 

We outline below the derivation-^ of the pseudospin 
hamiltonians and compatibility potentials in two spatial 
dimensions, for the square/rectangle, triangle/ /centred- 
rectangle, square/oblique, and triangle/oblique transi- 
tions, with number of variants Ny = 2,3,4 and 6 re- 
spectively. The 3D case is considered later. 



A. Square to Rectangle (SR.) hamiltonian: 

Nop =1,N v = 2 

Consider the prototypical square-to-rectangle or 'SR' 
transition, that is a two-dimensional analog of a tetrag- 
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onal to orthorhombic transition. For small distortions, 
the components of the symmetric Cartesian strain ten- 
sor are given by — 1/2 {d^u u + d^u^), where u(r) is 
the displacement vector and fj., v = x, y. We define linear 
combinations of the Cartesian components as three phys- 
ical strains, describing compressional (ei), deviatoric (e 2 ) 
and shear (63) distortions. 



- 2l i j. \ -21/ 

\&xx\&yy)i ^2 — 2 \~xx ' 



\ _ £3 , 

yy)i e3 ~ 2 ^ y ~ 



e yx)i 
(1) 



where c\ , C2 and C3 are symmetry-specific constants^. For 
the square reference lattice, C\ = c 2 = \[2 and C3 = 1. 
For the equilateral triangle reference lattice c\ — c 2 = 
C3 = 1. The pseudospin hamiltonian is obtained by the 
three steps given above, that we follow for all transitions. 
i) Scaled free energy, and compatibility potential : 
For the SR case, the deviatoric strain e 2 is the order 
parameter (OP). The compressional and shear strains are 
the non-OP strains. The scaled free energy is F = EqF 
where the overall Eq is an elastic energy per unit cell, 
and the dimensionless F = Fl + Fg + F non is a sum of 
three terms, 



Fnon — nAiei(r) — ^ —A, 



h(k)\ 2 . (7) 



r,i— 1,3 



k,i=l,3 



The scaled compressional and shear clastic constants can 
be expressed in terms of the (unsealed) elastic constants 
Cij in the Voigt notation, evaluated at To. For the cubic 
caseil,^i = (Cu + 2C 12 )/(Cu-C u ). The ratio A 1 /A 3 
is taken as fixed in simulations, for simplicity. 

For k = uniform contributions, the optimum non- 
OP strains are zero, at the parabolic fnon = minimum. 
For spatially varying k ^ contributions, the non-OP 
strains are to be minimized subject to the St Vcnant 
compatibility constraint ^ 1 1 > 1 2 that says distorted unit 
cells fit together in a smoothly compatible fashion, with- 
out defects like dislocations, so the displacement field is 
single-valued. The St. Venant conditions in the Carte- 
sian strain tensor e are n ' 12 (with ' T ' a transpose), 



V x (V x e) q 



0. 



(8) 



In two dimensions, the constraint in terms of physical 
F L = F G = J2 /c(Ve 2 ); F non = ^ /„ on ( ei , e 3 ). strains of Q is 



(2) 

where —¥ J d 2 r/a$ runs over all positions, and ag is 
a lattice scale for a computational grid. 

The dimensionless, scaled Landau free energy density 
in coordinate space is sixth order in e 2 (r) to give a first 
order transition^, 



-V 2 ei - -{dl - 9 2 )e 2 - -d x d v e 3 = 0. (9) 

C\ C2 J C3 



or in Fourier space 



O 1 e 1 (k) + O 2 e 2 (k)+O 3 e 3 (k)=0, 



(10) 



h{z2) = {T-l)el + e\{el-l) 
A scaled temperature is defined by 

t=(T- T C )/(T - T c ). 



(3) 



(4) 



There can be three 9/i/9e 2 = minima: at zero- 
strain e 2 = austenite, and at two martensite vari- 
ant minima of nonzero strain, e 2 = ±e(r). The order- 
parameter magnitude e at the variant minima is 



£ -(r) = [-(l + yT _ 3774)] 1 / 2 . 



(5) 



On cooling below the upper spinodal r = 4/3, two 
martensite variants appear ; they become degenerate 
with the austenite zero-state at r = 1 or T = To when 
e(t = 1) = 1; and for r < 1 the martensite wells become 
lower in energy. The austenite minimum disappears be- 
low the lower spinodal r = or T = T c . 

The cost of creating interfaces or domain walls is given 
by the usual Ginzburg term, with £ a wall thickness scale, 

Fg = ££ 2 (Ve 2 (rl) 2 = ]T £ 2 fc 2 |e 2 (fc)| 2 . (6) 



where the compatibility coefficients are 

— h 2 k 2 — k 2 Oh h 

Ox = —, 2 = HZ—OiL, Q 3 = x y 



Cl 



''2 



(11) 



The constrained minimization can be done through La- 
grange multipliers^, or by a direct substitution of the 
constrained solutionis e\ = — (0 2 e 2 + 3 e 3 )/0\ of (|10l) . 
into the non-OP free energy of (0, 

fnon = J2 \ [ A i K°2e 2 + 3 e 3 )/Oi| 2 + A 3 \e 3 \ 2 ] . 

(12) 

A free minimization in the remaining non-OP strain e 3 
determines it in terms of the OP e 2 . In fact, e% = 
-{O i 2 /A i )/[{pl/A 1 ) + {Ol/A 3 )} fori= 1,3. Substitut- 
ing into (|7|) yields the compatibility-induced interaction 
F C om P at{e2) = F non (ei,e 3 ), where 



Fcompat = ^ ^-U(k)\e 2 (k)\ 



(13) 



The compatibility kernel U in Fourier space is 



Finally, the non-OP strain energy is simply harmonic 
in compressional (ei) and shear (e 3 ) strains, 



AiUik) = u(k) 



0\ 



[{0\/A x ) + {Ol/A a )] 



(14) 
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or explicitly from (11), and c\ = c 2 = v2, C3 = 1, 

U{k) = u{k) \ k *~^ . (15) 

K < a A 3 K x K y 

Here the prefactor v(k) = 1 — 5g Q is inserted to make 

fnon ~ vanish for /c = uniform non-OP strains, as 
mentioned. In coordinate space, this is a powcrlaw in- 
teraction between OP strains U(R) ~ ^/R d , with sign- 
variation in angular directions yielding zero angular av- 
erage (52ftU(R) ~ f7(fc = 0) = 0), so it is not 'long- 
range' in the isotropic Coulomb ~ l/R d ~ 2 sense. The 
powerlaw anisotropic interactions are easily evaluated in 
Fourier space, and one need not resort to uncontrolled 
coordinate-space truncations to near-neighbor couplings, 
that may leave out some essential physics of the transi- 
tion. In coordinate space, 

F comp at(e 2 ) = \Ys a iW~ ?)e 2 (r)e 2 (7). (16) 
The formal partition function 

Z= /"n^2(r)exp(-/3F[e 2 (f)]), (17) 

T 

is dominated by free energy textural minima, that may 
be asymptotically found in a TDGL or rclaxational dy- 
namics, as done elsewhere 7 -, 

de 2 (r,t) _ 5F 

dt 5e 2 {r,t)' [ ' 

ii) Continuous strains to discrete pseudospins : 
One can approximate the partition function by retain- 
ing only the Landau-minima at fixed OP-magnitude val- 
ues \e 2 \ = e(r), and different OP signs (or in general, 
different angular directions of minima), while neglect- 
ing fluctuations about these minima. The continuous- 
variable strains are then replaced by discrete-variable 
pseudospins^ 

e 2 (f) -> e(r)S(r), (19) 

where the pseudospin has the three values S(f) = 0, ±1, 
to locate the minima at e 2 — 0, ±s(t). Although in zero 
stress the uniform austenite state is no longer a Landau 
minimum below the lower spinodal r = 0, the surround- 
ing nonuniform textures can exert local internal stresses 
to locally favor the zero value, even at low temperatures. 
Also, the free energy in OP strain always has a turning 
point at the origin to support dynamical transient ze- 
ros, that although few in number, could play a catalytic 
role in microstructural cvolutionAS. Hence we retain zero 
spin values at all temperatures, allowing their perma- 
nent/transient existence to be determined dynamically. 



With this substitution and S 6 = S 4 = S 2 = 1, 0, the 
approximated Landau free energy density at the minima 
can be written a.sH 

f L = e 2 (r)g L (T)S 2 (f); g L = (r - 1) + {f - l) 2 . (20) 

where e is in ([5]) . 

Hi) The reduced pseudospin hamiltonian : 
The partition function of (jTTJ) reduces to a sum over 
all the pseudospin configurations, with a temperature- 
dependent effective Hamiltonian in the Boltzmann 
weight, that can then be studied by the usual methods 
of statistical mechanics. Substituting (|T9"|) into the to- 
tal scaled free energy directly yields the hamiltonian in 
coordinate space, 

H(S) = F(e 2 -> eS), (21) 

where 

p H (S) = ^lJ2{9LS 2 (r) + e(VS) 2 } (22) 

r 

+ J2^-U(r-r<)S(?)Sm 

r.r' 

and £>o = 2£ , oe(r) 2 /T. This has the form of a gen- 
eralized spin-1 Blume-Capcl model 2 ^ as discussed later, 
but with temperature-dependent coefficients and power- 
law interactions. The hamiltonian is diagonal in Fourier 
spaced, 

PH=^Q (k)\S(k)\ 2 , (23) 

k 

where 

Q (fc) = D [g L (r) + £ 2 fc 2 + ^U(k)]. (24) 

B. Triangle/Centred Rectangle (TCR) 
hamiltonian: Nop = 2,Nv = 3 

Consider a two-dimensional crystal with equilateral tri- 
angles transforming to isosceles triangles, with three pos- 
sible such variants (Ny = 3), as there are three sides 
that can become the unequal side. The unit-cell changes 
from an equilateral triangle to a centred-rectangle. This 
'TCR' transition is the 2D version of the hexagonal to or- 
thorhombic transition observed in lead orthovanadate^. 
There are two order parameter a 4 ' 7 ' 17 ' 20 : the deviatoric 
strain e 2 , and the shear strain e^. The single non-OP 
variable is the bulk dilatation or compressional strain 
e\. Just like this TCR case, the square/oblique, tri- 
angle/oblique and cubic/tetragonal transitions also have 
the same OP (e2,e3) and a single non-OP strain ei, but 
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of course are distinguished by their different, transition- 
specific Landau polynomials, that induce different Ny 
directions of the vector pseudospins. 

i) Scaled free energy and compatibility potential : 
The free energy functional, invariant under the trian- 
gular point group symmetry, is 



F = £ (e2 ' es) + & (^2, Ve 3 ) + fnon (ei)}. 

r 

(25) 

The Landau free energy /j, for the TCR case describes 
the first-order phase transition between the single high- 
symmetry austcnite phase and the Ny = 3 martensite 
variants. It has a third-order term invariant under equi- 
lateral triangle symmetries, ^3 = e\ — 3e2e 2 . In scaled 
form, in coordinate space 



h(e 2 , e 3 ) = r{e\ + e 2 ) - 2(e| - 3e 2 e 2 ) + (e 2 + e 2 ) 2 . 

(26) 

Figure [T] shows the Landau free energy with three variant 
minima, for a low temperature. 



term in ((27)) vanishes. The three variant minima in the 
(e2,e3) plane form a triangle lying on a circle of radius 
s = e(r), where 



e(r) = -(l + ^l-8r/9). 



(28) 



On cooling below the upper spinodal r = 4/3, two 
martensite variants appear ; they become degenerate 
with the austenite zero-state at r = 1 or T = T when 
e(t = 1) = 1; and for r < 1 the martensite wells become 
lower in energy. The austcnite minimum disappears be- 
low the lower spinodal r = or T = T c . 

The Ginzburg term fa is quadratic in the OP-strain 
gradient so 

F G = £ € 2 (Ve,) 2 = J2 ek 2 \e t (k)\ 2 . (29) 



r,£=2,3 



k,£=2,3 



Finally, the non-OP term is simply harmonic in the 
single non-OP comprcssional strain, 

F non (ei) = = ^El^)! 2 - ( 3 °) 




FIG. 1. Contour plot in (e2,es) space of the Landau free en- 
ergy for the Triangle to Centred-Rectangle (TCR) transition, 
with parameters r = —2.5 and Eo = 1. The three degenerate 
energy minima correspond to the three martensite variants at 
this low temperature. 

In polar coordinates in OP space, the order parameter 
vector is e(r) = (e2,ea) = e(cos(f>, sin0) with magnitude 



The Landau free energy in polar 



e(r) ee (e| + e 2 ) 
coordinates, with 773 = cos 30 isi 7 - 

f L (e, 4>) = [(t - I)e 2 + e 2 (e - I) 2 ] + 2(1 - m )e 3 (27) 

The angular dependence is fr, ~ — cos 30. The mini- 
mum conditions df^/de = 0; df^/defr = yield four min- 
ima: at the e = austenite, and the three variant minima 
with sin 30 = where = 4> m = 0, 7r/3, 27r/3, so the last 



Substitution of the compatibility solution ei(fc) = 
— (C*2e2 + 03e3)/Oi as for (12) immediately yields the 
St Vcnant term F compat (e 2 , e 3 ) = F non (ei) in terms of 
the OP. 



F, 



compat 



E -y u w{k)e t (k)e e (ky 



k,l,l'=2,3 

The compatibility potential kernel 

U& = v{%)O l O v IO\, 
or explicitly from (fTT|) with c\ = Ci = C3 = 1, 



(k 2 ,-k 2 y ) 2 (2k x k y ) 2 
U22 — v — , U22 = v- 



(31) 



(32) 



k A 

2k x k y (k 2 — ky) 
U23 = v p — = U 32 . 



fc 4 



(33) 



In coordinate space, 
-F CO m P at(e2,e 3 ) = ^2 Uee>(r-r*)e i (r)e e ,(r'), 

f,r' £,£'=2,3 

(34) 

and as before, the powerlaw potentials fall off in 2D as 

U u > ~ i/R 2 - 

ii) Continuous strains to discrete pseudospins : 
For the TCR case (and other two-component OP 
cases), Ny > 3, but we do not simply get a generalized 
spin-j model with 2j + 1 states on a line, and j = Ny/2. 
Instead we obtain clock-like model s 17 ' 19 with discrete S 
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vector variables pointing to the polyhedronAV corners 
and centre, in A^p-dimensional space. Since the zero 
state is included, these may be termed 'clock-zero' Zjv^+i 
models^ 7 .. Note that, unlike pure clock Zjy models^ 7 -, the 
squared spin S(r) 2 = 1, is still a statistical variable and 
not a constant, because of the zero states. 

The continuous- variable strains at minima are replaced 
by discrete- valued pseudospinsi 7 ., 

© <») 

where the two components of the pseudospin have 
three variant values as in Fig [1] plus zero, 

a -(s)-(o)-(4)- (36) 

For the Ny = 3 variants, S = (cos0 m , sin</) m ) and S 2 
1 . with <j) m = 27r(m — l)/3 and m = 1,2,3. 

With this substitution and S 6 = S 4 = S 2 = 0, 1, the 
Landau polynomials again collapse into a simple form, 
bilinear in the pseudospins, 

Jl{t)=£ 2 9lS 2 ; g L { T ) = T-l + {e-\) 2 , (37) 

with £ as in (JSHJ) • 

Hi) The reduced pseudospin hamiltonian : 

In coordinate space the total pseudospin hamiltonian 

is 



PH = -f[ErE^3{9LSi(r) 2 +e\^S e (r)\ 2 } 

+ 4r £^=2,3 r*)^ W)L 

(38) 

and is a clock-zero modeli 7 - , with S having 3 + 1 

values of (f36|) . and with a compatibility kernel of ([33| . 
It is again diagonal in Fourier space, 

PH=^Y.Y.Q^'$)Se(k)S i ,(k)* 7 (39) 

k V 

with S(k)* = S(-k), as S(r) is real. Here 
Qo,w(k) = Do[{g L (T)+i; 2 k 2 }6t,e> + *±U U >$)\. (40) 

C. Square/Oblique (SO) hamiltonian: 

Nop = 2, N v = 4 

We consider the square/oblique or ' SO ' transition 
where the transition is driven independently by the devi- 
atoric e<i and shear e 3 order parameter strain a 7 ' 17 i 20 , as 
modified by a sufficiently strong coupling term. 

i) Scaled free energy, and compatibility potential : 



The Landau term has the scaled form 

h = r(e 2 + 4) - (4 - CJ2)(4 + ef) + 4(e« + e%) ~ C' 4 e 2 el 

(41) 

where C\ is a material-dependent clastic constant. In 
polar coordinates, with e = ( 62,63) = e(cos <f>, sin 0), it 

ft = [(t - l)e 2 + e 2 {e 2 - l) 2 ] + £ 4 (3e 2 - 2 + C'J2) cos 2 2^. 

. (42) 

The angular dependence is fi ~ cos 40. The five minima 
from dfL/de = O;3/t/90 = are the austenite zero 
state, and four variant minima with sin A<j) = in angular 
directions <fr = 4> m = 7r(2m— 1)/4 with m = 1, 2, 3, 4. The 
last term in (|42l) vanishes at minima, suppressing the C4 
material dependence. The four variant minima in the 
62) 63 plane for t < 4/3 form a square lying on a circle of 
radius e — e(r), where e is as in the SR case of 

ii) Continuous strains and discrete pseudospins : The 
strains at minima are replaced by pseudospins as in (|35[) . 
The discrete pseudo-spin has the five values^ 7 - 

For the Ny = 4 variant minima with S = 
(cos (/> m , sin <p m ) , and = TO = 7r(2m — l)/4 where 
m = 1, 2, 3, 4 , the spin magnitude is unity S 2 = 1. 
The Landau term becomes 

/l = £ 2 3l^ 2 ; 3 L = T-l + (e 2 -l) 2 (44) 

with e of (5). 

Hi) The reduced pseudospin hamiltonian: 
The Ginzburg and St Venant terms are the same as 
in the TCR case. The SO case clock-zero hamiltonian 
Z4 + i is formally the same as (|38|) . with S having 4 + 1 
spin directions of (|43[) , and the same TCR compatibility 
kernel of fl55J). 

D. Triangle/Oblique (TO) hamiltonian: 

Nop = 2, N v = 6 

The transition is, as in the TCR case, driven by a two- 
component OP2 i 17 i 20 e = (e2,e3). Here Ny = 6, so we 
need a square of the cubic term, I 2 to give six preferred 
angles. 

i) Scaled free energy, and compatibility potential : 
The scaled Landau free energy with up to sixth order 
invariants is^ 7 - 

h = (r - l)h + J 2 (/a - I) 2 + C 6 (/ 2 3 - / 3 2 ), (45) 

where I2 — e 2 = e 2 , I3 — e\ — 3e2e|, and C§ is a material 
constant. 
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In polar coordinates with 773 = cos 30, this ii 
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h = [(t - l)e 2 + e 2 (e 2 - l) 2 ] + C 6 e 6 (l - 77a 2 ). (46) 

The angular dependence is /l ~ — cos 60. Minimizing 
yields six martensitc variants with sin 60 = 0, at angles 
= 0m = 27r(m — l)/6 where m = 1,2, ..6, where the 
last term in (|46|l vanishes, suppressing the Cg material 
dependence. The six variants for r < 4/3 form a hexagon 
in the e2 , e3 plane, lying on a circle with radius s — e(r) 
of ©. 

n) Continuous strains to discrete pseudospins : 
With the usual approximation ([33)) of e(f) — > e(r)5(r), 
the pseudo-spin S(r) has seven values 



5 = 



±1 




1 V3 
- 1 - 2 



(47) 



The Landau term becomes 



J2 Ph MF {r) = J2 V(r)S(r) = £ V(k)S(k)*, (50) 



and C = h 



\ J2k Oo| (J (fc)| 2 - Here, V in Fourier and coordinate space 



2 
is 



and 



^ = [9l{t) + + ^U(k)\ a(k). (51) 



9dr)a{f) - e 2 V 2 a(rO + ^ E U (?~ r X r ") 



(52) 

The mcanhcld partition function is a product of local 
contributions 



f L (r) = e(r) 2 g L S(r) 2 ; g L (r) = (r - 1) + (e 2 - I) 2 , 

(48) 

with £ of ©. 

mj TTie reduced pseudospin hamiltonian 

The Ginzburg and St Venant terms are as in the TR 

case. The TO case clock-zero Zg+i hamiltonian is as 

in (|38[) with S having 6 + 1 spin values (|47|) . and with 

the compatibility kernel of (|33l) . 



III. LOCAL MEANFIELD IN TWO SPATIAL 
DIMENSIONS 



Z MF = £ e"^«- = LI E e-^ F « +c . (53) 

{S} f S(f) 

The self-consistency equation for the statistical average 
<x(r), with the constant C dropping out, is 



<j(r) = SWe-WWW J ^ ( 

S(f)=0,±l S(f)=0,±l 

that yields 



-0V(f)S(f) 



(54) 



With the pseudospin hamiltonians for SR, TCR, SO 
and TO transitions in hand, we now do local meanficld 
approximations^ for each of these cases. 



Square/Rectangle (SR) meanfield 

Nop = l,N v = 2 



cr(f) 



-2sinhF(r) 
1 + 2coshV(r)' 



(55) 



The equation can also be instructively obtained 
through the Gibbs-Bogoliubov inequality 



We write S(r) = a(r) + 5S(r), where cr(r) = (S(f)) 
is the spin statistical average, and substitute into the 
hamiltonian ([55)1 . Retaining only first order terms in 
SS(r) = S(r) — a(r), the meanfield hamiltonian is H = 
Hmf +0(SS 2 ). A similar approximation, with identical 
final results, can be done in Fourier space, with S(k) = 
a(k) + SS(k) substituted in ([39]) . 

The mean-field hamiltonian is then a sum of a local 
contribution and a constant, 



W. 



MF 



]0h MF (f)-C, 



(49) 



F < F var = Fo + (H — Ho)o, 



(56) 



where 



where the index refers to an average with a solvable ref- 
erence system iJo, taken here as H = — ^2pB(r)S(r). 
Here the local field B(r) is a variational parameter, and 
the free energy is F = -T^log[l + 2 cosh{/3 B(r)}]. 
The statistical average of S(r) in the reference system is 
a{f) = {S{-r)) a = 2sinh/3B/(l + 2 cosh ^B) and the av- 
erage of H — Ho can also be readily performed since the 
spins are uncorrelated. The optimal local field B(r), that 
minimizes F var through SF var /8B(f) = is then directly 
seen as B(r) = —V(r) with the same self-consistency 
equations as before. Hence V(r) is indeed the best molec- 
ular field to approximate the free energy of the original 
system. 
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The mean-field equations have been solved itera- 
tively under a cooling ramp in order to study long-lived 
glassy states^. Here we solve the equations for a fixed 
constant temperature r starting from an initial random 
configuration. With an input a(r) and an Fast Fourier 
Transform (FFT) to a Fourier a(k), it is easy to find V(k) 
from the definition (|5"Tj) . A reverse FFT to V(r) is used 
in (|55[) to obtain the next er(r), and the process repeats. 
Figure [2] shows twin microstructure obtained by solving 
the mean-field equations, with parameters as in the cap- 
tion. These twins are similar to those in experiment^, to 
relaxational simulations or to Monte Carlo simulations as 
in the inset. (Different phases, including certain maze- 
like textures are also seen in some parameter regimes^, 
but do not seem to appear in Monte Carlo simulations.) 

Thus a local mcanficld approximation to the pscu- 
dospin models is useful to study microstructure below 
ferroclastic transitions. 




2s(r) 



J2 h (r)S(r). 



(57) 



The mcanficld self-consistency equations become 



a(r) = 



-2sinh[y(r) - (D Q /2e)h(f)} 
1 + 2cosh[V(r) - (D /2s)h(r)] ' 



(58) 



Starting from random texture seeds with a small uniform 
external stress, h = 0.3, we obtain a uniform state of 
S = ±1 depending on the sign of h: the small stress 
picks out the global minimum. The twins are self-trapped 
metastable states, that are however quite rigid against 
stress: for a twinned initial state, a strong stress of about 
h = 4 needs to be applied to destroy the twins and to 
obtain the uniform ground state. Once the twins have 
vanished, the system fails to return to the original state, 
i.e. shows hysteretic behavior. 



2. Blume-Capel model phase diagram 

To make contact with treatments of the Blumc-Capcl 
model, we suppress the nonlocal couplings by setting 
A\ = 0. The Ginzburg term in ^ can be recast on 
a lattice, by setting the gradient to a discrete differ- 



ence operator V — > A, so that (VS")^ 



(Asy 



20 40 



00 120 



AS 2 -2Y, < ?r"> S{f)S(r'). Then the hamiltonian is pre- 
cisely a Blume- Capel model, bilinear in the spins (with- 
out the biquadratic term of the Blume-Emery-Griffiths 
model)22, 



FIG. 2. Microstructure obtained from the mean-field analysis 
of the Square-Rectangle (SR) spin model. The parameters 
are L = 128, £ 2 = 0.5, Eo = 3, and r = —2.5, and stiffness 
Ai = 4 with 2A1/A3 = 1. Twins are oriented along a diagonal 
as expected. Inset: Twinned ground state from Monte Carlo 
simulations for the same parameters. 



We now i) study effects of external uniform stress, ii) 
make contact with the phase diagram of the Blume-Capel 
model with uniform OP and iii) show how the Mean-field 
equations for cr{r) can be obtained through the least- 
action principle. 



1. Effects of external stress 

Twins with oriented, locked-in domain walls of posi- 
tive energy cost are metastable states, and the uniform 
single-variant state without domain walls is the global 
minimum in free energy. This can be seen by adding an 
external stress term h(f) with a simple linear coupling to 
the meanfield hamiltonian (|4T))) : 



A(r)5>(rf 



(59) 



There is temperature-dependence in the on-site crystal 
field term A(r) = A)(t)[ 5l (t) + 4£ 2 ]/(2/3), and in the 
ferromagnetic coupling J(r) = Dq(t)£ 2 / '(3. 

The model can be studied within the (uniform) mean- 
field approximation. An expansion of the mean-field free 
energy yields an analytical expression for the line of crit- 
ical points T c = zJ/3; and the location of the tricritical 
point, A c = |zJlog2, where z = 4 is the number of 
nearest neighbors. Figure [3] shows the well-known phase 
diagram of this model. Both A and J depend on the 
temperature r, so although a given temperature corre- 
sponds to a point, a cooling path is a line in the phase 
diagram. These lines intersect the first-order transition 
curve for r ~ 1, the Landau transition temperature be- 
tween ' paramagnetic ' austenite and ' ferromagnetic ' 
martensite. Figure 2] shows the meanfield phase diagram 
for two parameter planes {h,r) and (£, r). In a certain 
range of parameters, the spin model is consistent with 
the Landau theory that predicts a first-order phase tran- 
sition at t = 1. 
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FIG. 3. Mean-field phase diagram for the spin-1 model for 
the Square-Rectangle (SR) transition. The red solid line rep- 
resents second-order phase transition and the red dashed line 
represents a first-order phase transition^, meeting at a tricrit- 
ical point. The crystal field A(r) depends on the temperature 
r, so cooling is a phase-diagram trajectory. Three directed 
trajectories with different sets of parameters are shown, for 
cooling from r = 4/3 to r = —2. 5. The lines intersect the 
first-order transition line (dashed) for r ~ 1. 
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FIG. 4. Mean-field phase diagram of the Square-Rectangle 
(SR) transition in the (h,r) plane for £ = 1. The h = first- 
order phase transition temperature of r = 1 moves up with 
applied stress h. Inset: Phase diagram in the (£,t) plane. For 
£ large , there is a first-order phase transition line (dashed) 
with r = 1, while for £ small, there is a second-order transition 
line (solid), that moves to lower temperature. 



3. Field theory for o~(r) 



die point approximations of field theories. 

The partition function can be compactly written as 



Z = ^2 CX P ( o X! SrKrr'Sf" j : 
{S} \ fr' ) 



(60) 



where = /?J(r)c5 < ^ > - ^±U{r- f) - 2/3A{t)S^?, 
The first Kronecker symbol is non-zero only if r and 
r* are neighbors. We note that the kernel K??' can be 
recast using usual matrix notations 



Krr' = D (f 



e 2 V 2 - - 9L (r) 



(61) 



where (f\U\r*) = U(f — r 1 ). We may then 
use the standard Hubbard-Stratonovich transformation 

iB T A- l B tQ find 



dot A 



the exact integral representation of the partition func- 
tion 



Z = 



1 



(62) 



{2n) N / 2 VdetK 
with the action 

S[cp] = \Y, ct>(7)KZ r U(?) - J2 log (1 + 2 cosher))). 

fr' f 

(63) 

Finally, we define a(f) = J^-, K f },(j){r'). The partition 
function reads 



Z 



(64) 



where we have defined the formal measure 2?[cr(r)] = 
VdctK j (27r) jV / 2 n ? da (r). With this field-theoretical 
formulation of the partition function, our problem, the 
Mean-field approximation is obtained by minimizing the 
action 

S[v] = \Y i cr{r)K?riv{7) (65) 

ff' 

- J2 log ^ + 2 cosh(^ K^ai?))^ . 

Within the saddle-point approximation, o(r) = (S(r)) 
so that the field a is indeed the statistical average of the 
spin. We note that V(f) = -J2f' K PP ,(T (^) (<*!• <|52]))- 
so the minimization of the action yields the mean-field 
equation ([55]) . as expected. 



We show here how the partition function may be trans- 
formed to obtain a field theory for <j(f) = (S(r)), so the 
mean-field equation (|55|) results from a saddle-point ap- 
proximation of a functional integral. Other mean-field 
equations in this paper can similarly be obtained as sad- 



B. Triangle/Centred-Rectangle (TCR) meanfield: 

Nop =1,Nv = 3 

The TCR case spin hamiltonian is (|38|) . with spin val- 
ues S = (0,0),(l,0)(-l/2,±V3/2) of (J3BJ and g L as 



10 



in ((371). Since N OP = 2 for the TCR, SO, TO and 
CT transitions, their meanfield equations are all formally 
the same. From the substitution Se(f) = <f}(r) + SS(r) 
and linearization in SS(f) = Se(r) — di(r) the meanfield 
hamiltonian is (3H MF = /3/imf(^) — C, as in (49) but 
the local contribution is now 



For the TCR case sums over the Ny 
of d36j), this is 



02 



cosh(^T/ 3 ) 



2 [e~ v ' cosh(^) + C osh(^y 3 ) 



3 spin values 



(74) 



^/3W(r)= V e (r)S e (r) = Y. V ^)Se(k)* 



r,£=2,3 



(66) 

and c = §£(/3W) = |E^//( fc >*W , ■ Thc func - 

tions V 2 and V3 are defined in Fourier space by 

V 2 (k) = P^{k)o 2 {k) + ^-U 23 (k)a 3 (k), (67) 



V 3 (k) = Pg{k)a 3 (k) + £^U 32 (k)a 2 (k). 



where 



P°(fc) = D [g L + i 2 k 2 + ^-U u (k)}. 



(68) 



(69) 



03 



V3 



sinh(^V 3 ) 



2 [ e -% cosh(f ) + cosh( W 3 )] 



(75) 



where the position dependences of o~i (r) and (r) are left 
implicit. The coupled equations ([73)1. (|75|) were solved it- 
eratively on a i x L = 256 x 256 lattice with periodic 
boundary conditions with parameter values £ 2 = 0.8, 
A\ = 5, r = —6.5, and £b = 0.01. Here, and throughout 
the following other cases, To = 1.0, T c =0.9. 

Figure [5] shows the relaxed microstructure obtained af- 
ter 10 5 iteration steps. As in continuous- variable simula- 
tions in strains or displacements^ we also obtain nested- 
star patterns as observed in experiments^ for lead ortho- 
vanadate. However, unlike the continuous-variable mod- 
els which are computationally intensive, the spin models 
and the local meanfield solutions reach thc complex mi- 
crostructure relatively rapidly. 



Defining Pi in Fourier space as Pe(k) = Pg(k)ai(k), 
the coordinate space meanfield hamiltonian of (|66[) is 
then 



Square/Oblique (SO) meanfield: 

Nop = 2,N V = 4 



(70) 



1 £ U 23 {r- f)[a 2 {f)S 3 {f) + a 3 (?)S 2 (r))] 



The SO case spin Hamiltonian is formally the same 
as (H§ but with SO case spin values S = (±1/^2, ±1/2) 
of ([43]) , and is as in (|44|) . Doing a local meanfield ap- 
proximation as before, the formal self-consistency equa- 
tions of ([72]). ([TH]) become 



Thc partition function of this linearized meanfield Hamil- 
tonian can again be factorized as in ()53[) . 

The self-consistency equations as in ([5"3"[) . ([53)1 for the 
statistical averages {<?(r)} again have the constant C can- 
celling, so now with V — {V 2 , V 3 ) and S = (S 2 , S3), 



a(r) = J2 S(r)c 

S(r) 



V(r).S(f)^^2 c -V(r).S(r)^ 

S(r) 



(71) 



In terms of the Ny variant states 5* = (cos (f> m , sin <fi m ) 
with m = l,2...Ny this can be formally expressed for 
TCR, SO, TO and CT cases as 



02 



I + J2m-i e~( cos ^™ v ' 2+sin< ^ y3 ) 



(72) 



(T 2 



CT3 = 



2sinh(ii-^ 



V2 [1 + 2 cos 



S h(^)+2co S h(^)]' 



2sinh(i^; 



2 sinh( 



V2-V3 • 

V2 • 



V2[l 



2cosh(^; 



2 cosh( 



V2-V3 
V2 



)] 



(76) 



(77) 



The coupled equations were solved iteratively on a 
128 x 128 lattice with periodic boundary conditions and 
for different temperatures r, starting from an initial ran- 
dom texture. Figure [5] shows that the microstructure ob- 
tained for t = —2.5, has vortices, as in the classical clock 
models or in thc XY model. This vortex in the strain field 
differs of course, from an edge dislocation that is a struc- 
tural defect in thc displacement field. The pseudospin 
vortex at the meeting point of domain walls is character- 
ized by the winding number or topological charge 
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FIG. 6. Final state microstructure from meanfield equations 
for the Square- Oblique (SO). The four variants are in four 
different colors, and pseudospin spins orientations are also 
denoted. The discretized vortices (topological charge +1) and 
anti- vortices (topological charge —1), expected in the classical 
clock models, are also identified. Parameters used are L = 
128, £ 2 = 0.3, E = 0.2, and r = -2.5 and A x = 6. 



50 100 150 200 250 



FIG. 5. Final-state microstructure obtained from meanfield 
self-consistency equations for the Triangle-Centred Rectangle 
(TCR) transition: (a) strain component e2 ; (b) strain compo- 
nent e3. The color code of red (blue) corresponds to positive 
(negative) values, and green to zero. Note the sharp domain 
walls. Parameters are L = 256, £ 2 = 0.8, E = 0.01, scaled 
temperature r = —6.5 and stiffness A\ = 5. 

where 6(f) is the polar angle of the spin S(r), that equals 
4> m m the variant regions, and r, is an arbitrary con- 
tour surrounding the z-th vortex. The topological charge 
is qi = 1 for a vortex and g, = — 1 for an anti- vortex. 
Thanks to the periodic boundary conditions, we have 
J2i Qi = 0. Vortex solutions for complex fields at three- 
domain meeting points have been considered^ 3 -. 
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FIG. 7. Final state microstructure from meanfield equations 
for the Triangle-Oblique (TO). The three variants are in three 
different colors, and spin orientations are also shown. Only 
three of the six variants (see Fig 8), finally survive. Parame- 
ters used are L = 128, C 2 = 0.35, E = 0.2, scaled tempera- 
ture t = —2.6, and stiffness Ai = 5. 



Triangle/Oblique (TO) meanfield: 

Nop = 2,N V = 6 



D. 



The TO case hamiltonian is as in (|38|) but with TO 
case spin values S = (0, 0), (±1, 0), (±1/2, ±y/3/2) from 
(|47|). and gi, is as in (|48p. The general meanfield self- 
consistency equations of (|72l) . (|73| are then 



02 



2 sinh(F 2 ) + sinh(J) + sinh( J) 
1 + 2 cosh(V r 2 ) + 2 cosh(I) + 2 cosh(J) ' 



(79) 



0-3 



V3 



2sinh(7) - 2sinh(J) 



2 1 + 2 cosh(V2) + 2 cosh(I) + 2 cosh(J) : 



(80) 



where I = (V 2 + V3V 3 )/2 and J = (V 2 - V3V 3 )/2. Fig- 
ure[7]shows the ground state obtained from these coupled 
meanfield equations with parameters L = 128, £ 2 = 0.35, 
Eq = 0.2, t = —2.6 and A\ = 5. We note that discrete 
vortices at the junction of the six martensite variants, are 
seen only during the iterations through transient states 



12 



as in Figure [H The final state microstructure shows no 
vortices, and only three out of the six variants finally re- 
main, bounded by nonintersccting domain walls, as the 
other variants vanish during the course of the textural 
evolution. The suppression of vortices at least for these 
parameter values, could be due to the energy costs of the 
gradient and powerlaw terms. 
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FIG. 8. Transient state for the (TO) transition with the six 
variants with different colors, with parameters as for Fig 7. 



IV. PSEUDOSPIN HAMILTONIAN AND 
LOCAL MEANFIELD IN THREE SPATIAL 
DIMENSIONS 

We outline the hamiltonian derivations for the cu- 
bic/tetragonal case and then do a meanfield analysis. 
The approach can also be followed for other 3D tran- 
sitions^. 



A. Cubic/Tetragonal (CT) hamiltonian: 

Nop = 2,N v =3 

For the cubic-to-tetragonal or 'CT' transition, the 
symmetry-adapted strains are the dilatation e\ — 
(1 / \/S){e xx + e yy + e zz ), the two deviatoric OP strains 
e 2 = (1/V2)(e xx - e yy ), e 3 = (l/y/Q){e xx + e yy - 2e zz ), 
and the three shear strains = 2e yz , e§ = 2e zx , 

6Q 2iC X y. 

The OP components are the two deviatoric strains 
e = (63,62), and the remaining four non-OP comprcs- 
sional and shear strains are ei,e4,es,e6. The Landau 
free energy invariant under symmetries of the cubic unit- 
cell, was originally given by Barsch and Krumhansli^, 
where the cubic invariant is now I3 = (ejj — ie^e^), and 
in scaled form is 

h = t(4 + 4) - 2{e\ - 3e 3 e 2 ) + (e| + e 2 ) 2 . (81) 



The Ginzburg term is formally identical to (|29[) but in 
3D. 

The non-OP terms, harmonic in the four remaining 
physical strains are 

fnon = ^el + ^(el + el + el). (82) 

and are minimized subject to the compatibility con- 
straint ([5} in 3D. There are six equations, from cyclic 
permutations of the labels x, y, z of the two equations 

2d x d y e yz - d z e yy - d 2 e zz = 0, (83) 

d y d z e xx + d\e yz - d x d y e zx - d x d z e xy = 0. (84) 

By going to Fourier space one finds the second set is 
an identity, if the first set is satisfied. These con- 
straint equations can be recast in terms of the symmetry- 
adapted strains ei,e2...e6. Minimizing F non with these 
constraints (either through Lagrange multipliers^, or 
through direct solution for e4,es,e6 and minimization^ 
in the remaining ei), yields the non-OP strains in terms 
of the OP strains e 2 and e%. Substitution into the har- 
monic non-OP free energy yields the compatibility term 

F compa t = \ ^2^iU U '{k)e e (k)e e >(k)*, (85) 

i,i'=2,3 % 

where the kernels XJi y t> (k) in Fourier spaced are given in 
the Appendix. 

The procedure is formally just as in the TCR 
the CT case also has the same Nop = 2, Ny — 3, pseudo- 
spin values S = (0, 0), (1, 0), (-1/2, ±V3/2) and again 
e(r) = 3/4(1 + \Jl — 8r/9). The spatial dimension only 
enters in the 3D compatibility potential of kernels (|85|) , 
and in the 3D lattice positions r = (x, y, z) and Brillouin 
zone wave- vectors k = (k x ,k y ,k z ). 

B. Cubic/Tetragonal meanfield : Nop = 2, N v = 3 

We numerically solved the CT meanfield equa- 
tions (|74|) . (|75|) that are as for the TCR case, with the 
kernels as in the Appendix. We took a 32 x 32 x 32 lat- 
tice with periodic boundary conditions, and parameters 
C 2 = 10, E = 0.001, T = 1,T C = 0.9, and stiffnesses 
A\ = 4.8, A4 = 2.4. Fourier transforms enable a com- 
putation at each step, of the functions Vi2(r) and V^(r). 
Figurc[9]shows the microstructure, with twins at diagonal 
orientations, as found in continuous- variable simulations. 

V. OTHER RELATED MODELS 

Modified truncations of these structural-transition free 
energies can induce other hamiltonians, that can be stud- 
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FIG. 9. Twins in the (111) plane obtained from the meanfield self-consistency equations for the 3D cubic-to-tetragonal 
transformation. The color bar represents S v . The parameters are L = 32, £ 2 = 10, Eo = 0.001, scaled temperature r = —0.5 
and stiffnesses A\ = 4.8, A4 = 2.4. The microstructures (a),(b), and (c) show three different twin orientations obtained, for 
different runs. 



ied purely as interesting spin models in statistical me- 
chanics. 

Let us suppress the zero state, and fix only the cir- 
cle radius e — > e(t), while keeping all continuous polar 
angles, now denoted by 9(f), with values 2n > 9 > 0, 

so as in the XY model of planar spins, S = (cos 6, sin 9). 
The Ginzburg discrete-difference term of ^ then induces 
an XY- like ferromagnetic interaction. The Landau free 
energy in polar coordinates in all cases has angular de- 
pendence f L B cos N v <f> as in (J27J) , @2J), ([15]). 

Putting all this together, the free energy induces an XY 
ferromagnet model with a long-range potential, and a 
symmetry-breaking local field : 

PH = ^\-BY i <x»[N v 6{f)\ 

r 

-2e £ co S (6(r)-6(7))]+(3H c ({9}) (87) 

<r,r / > 

Here as S 2 = 1, there is no quadratic local term, and 
/3Hc({0}) is a term coupling the continuous-angle vari- 
ables cos 9(f) and sin 0(f), 

PHc[{9}\ = U ™ cos9 ^ cose (?) 

r,r' 

+ C/33 sin 9(f) sin 9(f) 

+ U23 cos 9(f) sin 9(f) + U 32 sin 9(f) cos 9(f) 

(88) 

and a partition function 

Z= [ ld9(f)cxp(-(3H{{9(f)};T}). (89) 

J [0,2^ p 



This model includes all angles, even away from minima, 
and so can describe slowly transiting states across Ny 
saddle points, as in experiment 2 ^. 

Similar XY models with symmetry-breaking fields 
(without the powerlaw interaction) have been studied 2 ^. 
A dual transform in that case extracts the topological 
vortices with logarithmic interactions, and in this model 
could also induce a powerlaw anisotropic vortex interac- 
tion. A real-space renormalization group analysis of the 
2D Coulomb gas as in the Kosterlitz -Thouless transition 
is well-known 2 ^, and could be repeated for this model. 
Renormalization flows in the context of martensitic tran- 
sitions have been studied in other models 2 ^. 

For strong symmetry-breaking in minima angular di- 
rections ( \B\ large), the continuous angle 9 become dis- 
crete and takes on values that we denote as 6 <f> = 4> m , 
and one gets pure clock models (Zjvv) with a Hamilto- 
nian that now has powerlaw potentials, 

maock = ~D o e cos[0(r)-0(r')]+/3i/ c ({0}). 

<rf> 

(90) 

We can even make one more approximation by reducing 
the XY interaction to a Kronecker-delta coupling, yield- 
ing a g-state PottsiS model with q = Ny. 

<rr'> 

Potts hamiltonians with large number of spin compo- 
nents q have been studied as models for configurational 
glasses^. 

VI. CONCLUSION 

A standard approach to obtaining microstructure of 
structural transitions is to solve evolution equations for 
relaxation to a minimum, in continuous variables such 
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as displacements, phase fields or strains^. We have 
here considered the reduced hamiltonian models in dis- 
crete pseudospins describing four structural transitions 
in two dimensions, as well as the three dimensional 
cubic-to-tetragonal transition. These 'clock-zero' mod- 
els have a zero state as well as clock states, and the 
pscudospin hamiltonian has an on-site term, an exchange 
interaction and a powcrlaw interaction term. For the 
square/rectangle case, the pseudo-spin model without 
powerlaw interactions corresponds to the Blume-Capel 
spin-1 model with temperature-dependent couplings. Us- 
ing a local meanfield approach, we have obtained the mi- 
crostructure for 2D and 3D transitions, as obtained in 
continuous-variable strain dynamics. For example, the 
characteristic nested star microstructure of the triangle 
transition emerges easily from the meanfield solution. 

The textures of the squarc/obliquc (SO) and trian- 
gle/oblique (TO) transitions, with Ny = 4, 6, which 
have not been previously studied, include vortex con- 
figurations of the Zjv v _|_i clock models, at intersections 
between variant domain walls. The SO final microstruc- 
ture has positive/negative vortices in regular patterns, 
and all four variants are present. For the TO case, at 
least for particular parameters, we find the six-variant- 
vortices appear only as transient solutions, with the fi- 
nal state having no vortices, with only non-intersecting 
closed-domains of three variants. Finally, for the three 
dimensional cubic/tetragonal transition, we obtain the 
diagonal twinning that is consistent with previous stud- 
ie a 4 i9'7 . In all cases, the local meanfield final microstruc- 
ture emerges relatively rapidly, compared to the slow evo- 



lution towards steady-state of the continuum differential- 
equation dynamics. 

Further work can involve studies of pscudospin hamil- 
toniansii for other structural transitions in 2D and 3D 
in the local meanfield approach. By including quenched 
disorder, such pseudo-spin models may be used to study 
strain glass behavior in martensitic alloys^, and relate 
solutions to the tweed precursors^ in analogy with spin- 
glass like behavior, and to random-field models^. Monte 
Carlo simulations can be used to study martensitic nuclc- 
ation and growth^. Other related spin models of interest 
in their own right may include geometric nonlinearitics 
that yield complex heirarchical-twin patterns^. 

In conclusion, the discrete-variable pscudospin model 
in local meanfield approximation, is therefore a useful 
approach to the study of martensitic texturing. 
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tetragonal transition. To do so we define the coefficients 

( s) 

Oa and O a by 



(Al) 



i(5) _ _1 n„2 , n(5) _ - I n (5) _ 1 i 



Or = -t={< + K),0\ 



V3 



V2 



y0 3 — /^(^k x k z ), 



V6 



(A2) 



0[ 



v/2 



(k x ~k y ),0 3 ——j=(k x + k y ), 



O4 — kyk z ,0^ — k x k z , Og — k x k. 



(A3) 



(A4) 



Let (5i s) = { a s) /O s and G a p = ^(As/A^O^ . 
The compatibility kernel for the cubic/tetragonal transi- 
tion can then be written as the 2x2 matrix 



Appendix A: Kernels for the cubic-to-tetragonal 
transition 

In this Appendix we state the explicit form of the bulk 
kernels Uu> obtained elsewhere^ for the 3D cubic-to- 



Utr = v(k) — — . (A5) 

1 + Crll 



where v(k) = (1 — 6% ) sets the non-OP harmonic-energy 
contribution for uniform strains to its minimum value of 
zero. 



